Pauli exclusion operator and binding energy of nuclear matter 
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Brueckner-Hartree-Fock calculations are performed for nuclear matter with an exact treatment 
of the Pauli exclusion operator in the Bethe-Goldstone equation. The differences in the calcu- 
lated binding energy, compared to the angle-average approximation, which is commonly used, are 
non-negligible. These difference exhibits a specific density dependence, which shifts the calculated 
saturation point towards smaller densities. This effect is observed for various versions of modern 
models for the NN interaction. 



It is one of the very central and very old projects of nuclear structure theory, to evaluate the saturation properties 
of nuclear matter from a realistic nucleon-nucleon (NN) interaction without any adjustment of a free parameter. The 
^ ■ so-called hole-line expansion or Brueckner-Bethe-Goldstone theory has been one of the tools for solving this many- 
body problem, which has already been used for many years Jl]-§]]. These early investigations were successful to some 
extent. The inclusion of NN correlations in the lowest order of the hole-line expansion, the Brueckner-Hartree-Fock 
t^. ' (BHF) approach, turned out to be very important. Realistic models of the NN interaction like the Reid soft-core 
potential |^| yield an energy of nuclear matter around 150 MeV per nucleon if the effects of correlations are ignored 
in a mean-field or Hartree-Fock calculation. The BHF approach provides a drastic improvement leading to an energy 
per nucleon of —11 MeV, which was only by 5 MeV off from the empirical value of —16 MeV per nucleon. Attempts 
have been made to improve the description of the saturation point further by exploring different NN interactions. 
It turned out, however, that BHF calculations using these various NN interactions yield results for the saturation 
^vq ■ point, which fall on the so-called Coester band ||. They either predict too small binding energy at the empirical 
y—\ \ value for the density, or about the correct energy at a density, which is too large by a factor of two, or results in 
OO > between. Comparison of BHF with variational calculations furthermore demonstrated that the inclusion of three-hole 
line contributions seems to be necessary to obtain a reliable estimate for the binding energy of nuclear matter [0J|]. 

During the last years some progress has been made in this field. It has been shown that a continuous choice for the 
particle spectrum |],|l(j (see also discussion below) accounts for the main part of the effects of three-body correlations. 
The discrepancy between the calculated saturation points of nuclear matter and the empirical one has significa ntly 
been reduced by considering relativistic effects within the Dirac-Brueckner-Hartree-Fock (DBHF) approach [[ll] |l3| . 
Finally, it should be mentioned that a new generation of realistic NN potentials has been developed [p^4l6|, which 
yield very accurate fits of proton-neutron and proton-proton scattering. These new potentials, which are essentially 
phase-shift equivalent remove a large part of the discrepancies observed between older models of the NN interaction 

;_i ' Because of these improvements, the time seems to be appropriate to test the reliability of approximations which 
are generally employed in the BHF calculations of nuclear matter. One of the central equations to be solved in the 
BHF approximation is the Bethe-Goldstone equation, which we may write in the momentum representation 

G s ,T(k,k',K,u) = Vs, T (k,k') + f d 3 P Vs, T (k,p} Q{Pl * ] G s ,t(p, k', K, u) (1) 
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In this equation K represents the center of mass momentum of the interacting pair of nucleons while k, k' and p stand 
for relative momenta, which are related to the single-particle momenta according to 

K=^{k 1 +k 2 ) and fc = i(*2-*i), (2) 

uj denotes the starting energy, while Vs,t an d Gs.t refer to the matrix elements of the bare interaction and G-matrix 
respectively for two nucleons with total spin S and isospin T. The operator Ho is used to define the energy spectrum 
of the intermediate two-particle state (K,p). The conventional choice for these energies of particle states above the 
Fermi surface has been to replace Hq by the kinetic energy of these states. In the so-called continuous choice, one 
assumes that the single-particle energies for these particle states as well as for the hole states are calculated from the 
kinetic energy plus a mean field contribution, which is calculated in a self-consistent way from the G-matrix by 
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The second line exhibits a parametrisation of these single-particle energies in terms of an effective mass and a constant 
potential, which is often used. The Pauli operator in the Bethe-Goldstone equation ([!]) prevents scattering into 
intermediate states with momenta p\ = K — p and p% = K + p, which are smaller than the Fermi momentum Uf- 
Therefore the value of Q(p, K) depends on the angle £1 between the center of mass momentum K and the relative 
momentum p. Employing a partial wave expansion of the two-particle states, the matrix elements for this Pauli 
operator can be written ]18| 
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with O(x) defining the step function. One finds that the Pauli operator is diagonal with respect to the modulus of the 
momenta K and p, the spin S and the projection quantum number m. It has non vanishing matrix elements between 
states of different I and J with the restriction of parity conservation (I + 1' must be even) and its value depends on the 
projection quantum number M. This implies that also the solution of the Bethe-Goldstone equation is not diagonal 
with respect to the angular momentum quantum numbers I and J and depends on M, 

(klJ\G ST M(u,K)\k'l'J'). (5) 

In order to simplify the calculation this Pauli operator is usually replaced by the so-called angle-average approximation 
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The angle-average approximation yields matrix elements for the Pauli operator which are diagonal in the I and J and 
independent on M. This means that the Bethe-Goldstone eq. (|l|) can be solved separately for each partial wave and 
the resulting G-matrix will be diagonal in J and independent on M . 

The matrix elements of G can then be used to evaluate the total energy per nucleon. In the case of the exact Pauli 
operator this energy is given as 
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with f2 the angle between the direction of the relative momentum k and the center of mass momentum K. In the 
calculations discussed below we consider the coupling of partial waves up to J, J' < 6. For angular momenta larger 
than 6, the Born approximation is used. If the matrix elements of G are diagonal in the total angular momentum 
( J = J') and independent on the projection quantum number M, as it is the case in the angle-average approximation 
for the Pauli operator, this expression can be rewritten into 
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which corresponds to the standard expression discussed e.g. in 0. 

Results for the calculated binding energy per nucleon as a function of the Fermi momentum are displayed in Fig. [j], 
considering various NN potentials. The continuous choice (figures in the upper half) as well as the conventional 
choice were used for the single-particle spectrum. For each density, NN interaction and particle state spectrum, a self 
consistent single-particle spectrum has been determined using the effective mass parametrisation defined in eq. (|^). 
The total binding energy has then been calculated using the angle-average approximation for the Pauli operator 
(dashed lines) as well as the exact treatment (solid lines). 

A general feature can be observed, which is independent on the NN interaction and the choice for the single-particle 
spectrum: The angle-average approximation tends to underestimate the binding energy per nucleon at low densities 
but overestimates it at higher densities. The effects of the exact treatment of the Pauli operator is not very large 
around and below the empirical value for the saturation density. This is in agreement with older studies of the 



angle-averaged approximation [TL8 20 1, in which matrix elements of G were compared at those small densities. The 



characteristic density dependence for the exact treatment of the Pauli operator, however, leads to a non-negligible 
shift in the calculated saturation, i.e. the minimum in the energy versus density curve. This correction moves the 
calculated saturation points to smaller densities and smaller energies. It is worth noting that the saturation points 
calculated for the Bonn C potential (E= —14.97 MeV at fcp=1.42 fm^ 1 ) and the neutron-proton part of Argonne Vis 
potential (E= —13.89 MeV at fci?=l-37 fm _1 ) are rather close to the empirical value. 

The effects of an exact treatment of the Pauli operator are larger for the continuous choice of the single-particle 
spectrum than for the conventional one. This is quite plausible as the continuous choice yields a larger sensitivity to 
the proper treatment of states around the Fermi momentum. A similar argument can be used to explain the fact that 
the Pauli effects are larger for NN potentials which have a slightly stronger tensor force (Bonn C and Argonne Vig) 
than is observed for the Bonn B potential, which contains a weaker tensor component. 

In summary we would like to point out that the effects of an exact treatment of the Pauli operator in the Bethe- 
Goldstone equation are not very dramatic in particular for nuclear matter at small densities. A characteristic density 
dependence of these Pauli corrections, however, lead to a non-negligible improvement in the calculated saturation 
points. Therefore many-body calculations going beyond the BHF approximation should take these effects into account. 
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FIG. 1. Calculated binding energies per nucleon for nuclear matter as a function of the Fermi momentum fcf - Results are 
presented for the angle-average approximation of the Pauli operator (dashed lines) and the exact treatment (solid lines). The 
potentials Bonn B and C, defined in |l9| , as well as the neutron-proton part of the Argonne Vi8 potential were used for the 
NN interaction. The continuous choice for the particle state spectrum in the Bethe-Goldstone eq.(jjj) has been used to obtain 
the results displayed in the upper part, while the conventional choice has been used to calculate the results shown below. 
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